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0^ , The spectrum of the evolution operator associated with a nonlinear stochastic flow with additive 

0^ ■ noise is evaluated by diagonalization in a polynomial basis. The method works for arbitrary noise 

strength. In the weak noise limit we formulate a new perturbative expansion for the spectrum of 
O , the stochastic evolution operator in terms of expansions around the classical periodic orbits. The 

D ' diagonalization of such operators is easier to implement than the standard Feynman diagram pertur- 

|Jh \ bation theory. The result is a stochastic analog of the Gutzwiller semiclassical spectral determinant 

p^ ■ with the "?i" corrections computed to at least two orders more than what has so far been attainable 

f^ ' in stochastic and quantum-mechanical applications, supplemented by the estimate for the late terms 

in the asymptotic saddlepoint expansions. 
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I. INTRODUCTION 



The periodic orbit theory relates the spectrum of the Fokker-Plank operator and its weighted evolution operator 
generalizations to the periodic orbits via trace formulas, dynamical zeta functions and spectral determinants P,H. 
_ For quantum mechanics the periodic orbit theory is exact on the semiclassical level [ [L3[ , whereas the quintessentially 
.*^ quantum effects such as creeping, tunneling and diffraction have to be included as corrections. In particular, the 
-Y-\ higher order h corrections can be computed perturbatively by means of Feynman diagrammatic expansions p| . 

I . The evolution operator formalism allows us to calculate long time averages in a chaotic system in terms of the 

p\J ' eigenvalues of evolution operators. The simplest example is provided by the Perron-Frobenius operator 
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Cp{x') = / dx5{f{x)-x')p{x) 



for a deterministic map f{x) which maps a density distribution p{x) forward in time. Our purpose here is to develop 
effective methods for computation of spectra of stochastic evolution operators. In case at hand, already a discrete 
time 1-dimensional discrete Langevin equation |0,|[, 



t-n+l 



/(x„) + a£,n , (1) 



^ ' with ^„ independent normalized random variables, suffices to reveal the structure of the dependence on the noise. 
. . 1 We treat a chaotic system with external noise by replacing the deterministic evoluton (5-function kernel by the 

Fokker-Planck kernel corresponding to (0) , a sharply peaked noise distribution function 

C{x\x)^6M{^)-^')- (2) 

In ref . |l| we have treated the problem of computing the spectrum of this operator by standard field-theoretic Feynman 
diagram expansions; in ref. [gj we offered a more elegant formulation of the perturbation theory in terms of smooth 
conjugacies. This time we evaluate the evolution operator in an explicit polynomial basis. The procedure, which 
is relatively easy to automatize, enables us to go three orders further in the perturbation theory; in the language 
of Feynman diagrams of ref. ||l|, in the new approach developed below we are able to compute the perturbative 
corrections up 5-loop level, as well as study the asymptotics of late terms in the perturbative expansion. 

The paper is organized as follows: in refsectFLOWS we review the evolution operator formalism for smooth flows. 
In sect, p wc explain the theorems that guarantee that spectral determinants for Axiom A systems are entire. In 
refsectNUMERICAL we discuss the numerical tests of our perturbative expansions. 



II. MATRIX REPRESENTATION OF Perron-Frobenius operator 

We shall sketch here the basic ideas behind the proofs that the classical spectral determinants are entire, without 



burdening the reader with too many technical details (rigorous treatment is given in refs. |19 20|JlO(| ). The main point 
is that the spectral determinants are entire functions in any dimension, provided that 

1. the evolution operator is multiplicative along the flow, 

2. the symbolic dynamics is a finite subshift, 

3. all cycle eigenvalues are hyperbolic (sufficiently bounded away from 1), 

4. the map (or the flow) is real analytic, ie. it has a piecewise analytic continuation to a complex extension of the 
phase space. 

As in physical applications one studies smooth dynamical observables, we restrict the space that £ acts on to smooth 
functions. In practice "real analytic" means that all expansions are polynomial expansions. In order to illustrate how 
this works in practice, we first work out a simple example. 

A. Expanding maps with a single fixed point 

We start with the trivial example of a repeller with only one expanding linear branch ^ 

fix) ^Ax |A| > 1 . 
The action of the associated deterministic, noiseless Perron-Frobenius operator is 

£0(y) = J dx5{y - Ax)0(x) = ^'/'(zz/A) . 
From this one immediately identifies the eigenfunctions and eigenvalues: 

^^"^^^"' - = 0,1,2,... (3) 

We note that the eigenvalues A~"^^ fall off exponentially with n, and that the trace of C is given by 

1 °° 1 1 

*'^^ MS"^"" ^ |A|(1-A-^) ^W^V ^^^ 

A similar result is easily obtained for powers of C^ and for the spectral determinant one obtains: 

det(l-z£) = n(l-^)=EO^-*'' 

t = — z/|A|, where the cumulants Qk are given explicitly by the Eulcr formula [M 
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(5) 



These coefficients decay asymptotically faster than exponentially, as A"'^'^'^"^^/^. An intuitive way of comprehending 
the spectrum is to view the composition with an expanding map as a "smoothing" operator; fast variations in the 
initial density distribution (/)(a;), corresponding to high powers of a:, are wiped out quickly by the smoothing operator. 
In a suitable polynomial basis 4>n{x) the operator has an explicit matrix representation 

OC 

(£0)„(x) = ^ Lnm<i)m{x) . 
m=0 

In the single fixed-point example (H), 0„ = y", and C is diagonal, Cnm = '5m„A~"/|A|. In general case, a matrix 
representation can be constructed by means of Cauchy complex contour integrals. 

The simplest example of how the Cauchy formula is employed is provided by a nonlinear inverse map ip — /~^, 
s — sign('0') 



C(f){w) — dx S{w ~ f{x))(f>{x) — s il}'{w) 4>{tp{w)) . 



Assume that i/) is a contraction of the unit disk, i.e. 

|V'(w)| < 9 <1 and W {w)\ < C < oo for |w| < 1 , 
and expand in a polynomial basis by means of the Cauchy formula 

dw 4>{w) 



4>{x) = ^ X"(j)r, 
n>0 



27ri w — X 



In this basis, £ is a represented by the matrix 

£(l>{w) = Y, w"i 

Taking the trace and summing we get: 



mnYn ; 



dw 4){w) 
27ri «;"+-'- 



dw s 'ip'{w){'ip{w)Y' 
27ri 



w 



,?n+l 



(6) 



(7) 



tvC = J2Lnn= fl^ 



n>0 



dw S 'lp'{w) 



2711 w — ^j(w) 



This integral has but one simple pole at the unique fix point w* = ^j{w*) ~ f{w*). Hence 

siP'{w*) 1 



tiC = 



i-^'{w*) \nw*)-i\ 



in agreement with ^. 

In practice we do not evaluate the matrix elements (M) by Cauchy integrals; instead, we use the polynomial ba- 
sis 1 10 1^, with the left basis given by derivatives t^^tt, and the right basis by polynomials ■§-^- An analytic function 

can be written as F{y) = ^Fm^, and the coeficients Fm can be obtained by F„ 



d"^F{y) 



dy" 



or in other words, 



y=o 



the dual basis of ^ is ^. So the Lik matrix element is obtained by acting the operator on ^ and differentiating 
the result I times at y = 0: 
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Qylm 
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CsAy'^y) 

5{y' + x^+l - f{y + x^)) 



V 



(8) 



While it is not at all obvious that what is true for a single fixed point should also apply to a Cantor set of periodic 
points, the same asymptotic decay of expansion coefficients is obtained when several expanding branches are involved. 
For a two-branch repeller, the procedure is the same, except the integral (0) picks up a contribution from each branch. 

From bounds on the elements Cmn one verifies |l^,g^,|l^ that they again fall off as A"*^ /^, concluding that the C 
eigenvalues fall off exponentially for a general Axiom A one-dimensional map. 



III. Evolution operator COMPOSITION FOR PIECEWISE- ANALYTIC EXPANDING MAPS 

Suppose we can decompose the evolution operator in a sum of two operators C — d^ + Ci. Here we shall study a 
noisy map with a binary Markov partition M = {Alo, A^i}- We distinguish two branches of the map : /(x) = fo{x) 
if X e Aia, and f{x) = fi{x) if x G Ali. In the case of Gaussian noise the corresponding operators are 



Co{x',x) = 



1 



2na 



.f.-^i^'~fo(^)f 



£^(^',^)^^e-i^(^'-/^(^»' 



2Tra 



X e Mo 
X e Ml . 



(9) 



If you visualize C as an [2n x 2n] matrix, in the weak noise limit the matrix elements for whom x is in the other 
partition are negligeable and can be set to zero, so Lq and £i are [n x 2n] matrices. 

For piecewise- analytic maps the decomposition is exact; for stochastic operators is assumes that the overlap of the 
two kernels is insignificant and can be neglected. Then we can write: 

lndct(l - z{Cq + £1)) = tr ln(l - z{Ca + £1)) 

= -V-tr(£o+/:ir. (10) 

^ — ^ n 



n 

Tl=l 



A contribution to this sum corresponds either to a prime cycle p = S1S2 • • • s„ , Si £ {0, 1}, or to its repeat 

tr -Cpr tr i^i-'S^ ^S2 ' ' ' '^Srip ) 

Expanding the 71th power we get contributions from all symbol sequences of length n: 

00 
tr (£0 + CiT = II ^P II '5„,„^.tr {CpY, (11) 

p r—l 

where p denotes a prime cycle itinerary composed of O's and I's, Up is the prime cycle length and r is the repetition 
number. For instance, p = Oil then £011 = £o'Ci£i = CqC^. The cyclic property of trace yields the Up factor. 
Substituting into (llO) we obtain 



lndet(l-z(£o + £i)) = -mi trC;. (12) 

p r—l 



The sum over repeats yields a factorized formula of the dynamical zeta function type 

det(l - z{Cq + Ci))=Y[ det(l - 2""/:^) . (13) 

p 

The operators are defined on piecewise monotonic maps, so there is only one periodic orbit for a given prime cycle 
itinerary. 

IV. SADDLE POINT EXPANSIONS IN TERMS OF PRIME CYCLES 

In the weak noise limit the kernel is sharply peaked, so it makes sense to expand it in terms of the Dirac delta 
function and its derivatives: 






2 ^3 



where 



= 5{y) + a2y<5(2)(y) + a^^S'^^^y) + .... (14) 

5^^\y)^^5{y), 

and the coefficients «„ depend on the choice of the kernel. We have omitted the ^'■^•' {y) term in the above because in 
our applications we shall impose the saddle-point condition, that is, we shift / by a constant to ensure that the noise 
peak corresponds to y = 0, so 5^{Q) ~ 0. For example, if 5„{y) is a Gaussian kernel, it can be expanded as 

= %) + Y'5(2)(y) + ^j(4)(y) + .... (15) 



V. Evolution operators IN A MATRIX REPRESENTATION 

In this section our goal is to calculate noise corrections to the leading eigenvalue of the Perron-Frobenius operator. 
Expression ( |lO| ) showes, that in order to do that first we have to calculate tiC", and from equation ( [ll| ) we see that 
tr £" should be generated from traces of C on periodic orbits. This section is a breaf review of how this was carried 
out in practise using the matrix representation of the Perron-Frobenius operator. 

If the coordinates of a prime cycle are xi, ..., a;„ , the operator of a periodic orbit segment is 



^sAy\y) = Yl 



-<5™(y' + x,+i-/(x.+y)). 



(16) 



m=0 



where Xn +i = xi and am are the moments of the noise. The full contribution to the trace from this periodic orbit 



is: 



tr; 



The calculation of this contribution can be computerized if we represent (Hq) in a matrix form using the polynomial 
basis (^). The matrix elements of Cs^ are: 



QT. 



dy' 



CsAv'^y) 



y" 



E 



anO_^(i) 



I m-^-n^m' ' 



/ n^raax{m' —ra.Q) 

where the B matrix is the representation of the noiseless operator: 



^kk' 



Qk 



dy" 



5{y' + x^+i - f{y + x^)) 



yk^ 

k'\ 



(17) 



(18) 



If the Dirac-delta in B acts on "^ we get: 



, ., , ^ fl ^ ^^/■' {! ^{x^ + l+y')~Xi) 

dyo[y +x,+i- f{xi+y))— - 



k" 



sign(/0 d 
{k' + iy.dy 



-{f-\x,+i+y')-x,) 



k'\ k'\\nj-^{x,+,+y'))\ 

k' + l 



(19) 



where sign(/') is the sign of /'(/^^(x^+i + y'j). There is no summation for the branches of the prepimages of /, since 
each branche takes us back to a differernt surrounding, so at each point the branch is choosen by the previous point 
of the orbit. The k, k' matrix element of B is the kth derivative of this at y' = 0: 



Bl 



sign(/') d- 



,fc+i 



'=''=' (fc' + l)! dy"'+ 
F\y') = r\x,+,+y')-x. 



iU/\k'+l 



TF\y' 



(20) 



y'=o 



Let us write F^{y') in power series of y' . li k' > k then after the derivates are taken every term in Bk^k' has at least 
a {y')'^ ~^ coefficient. We have to evaluate the derivative at y' = 0, so we see that i? is a triangular matrix, B^.k' = 
if k'>k. 



A. Multinomials 



The non-zero matrix elements are p3| 



a=i 



(21) 



where the sum (^) goes over all non- negative integers such that: 

fli + 2a2 + ... + nan = n , ai + a2 + ... + an ^ m (22) 

and the multinomial coefficient is: 

n\ 

(n\aia2...an) = -r^s tttt: ; — -, — r^ r (23) 

^ ' ' (l!)"iai!(2!)°2a2!...(n!)'»"a„! ^ ' 

If we expand F^ijj') in a Taylor series, the constant term is zero because of f~^{xij^i) — Xi. So we can write: 

^'(y') = ^ + f:f/ (24) 

1=2 

We apply the formula ( |2l| ) to F^{y') with power fc' + 1: 

{F^{y')f+'^{k' + iy. Y. y—Y^[n\a,a2...an)'j^^{F^2r...{F:r-. (25) 

n=k' + l 

In the (fc + l)th derivative of this aX y' — Q only the n = fc + 1 term is non-zero. So considering only one branch, with 
positve sign of /^: 

Blw = E(^ + lWia2-au+iy^^{FiY'-{Fl+iT''^' (26) 

ai + 2a2 + ...(fc + l)afe+i =fc + l (27) 

ai+a2 + ...afe+i =fc' + l (28) 

(fc + l|a,a2...a.^.,)' = (i,).,,^,(2,).,j^t!(('fc^+ !)!)-..«,+,! ^'^^ 

(Here F^ is the Zth derivative of F' on the particular branch which we have to choose.) For the diagonal and the 
nearest off-diagonals the matrix elements read as: 

Bm,n = ^,n+l ' "1 = 0,1,... (30) 

_ (m + 2)(m+l) f,' 

^™+i.™ 2 A™ ^ ' 

_ (m + 3)! {F^f , (m + 3)! fj 

(m + 4)! (fj)3 (m-^4)! f^fj (m + 4)! Fi 

"+^'" 48(m-2)!A™-2 ^ 12(m-l)!A™-i 24m! A™ ^ ^ 

etc... , 



B. Noiseless case 

In the noiseless case the B' matrices are the representations of the Perron-Frobenius operator. Since they are 
triangular, their eigenvalues are: 

Am = Bmni = j:;;^ ■ (34) 

Multiplied, the triangular matrices yield a triangular matrix, and the diagonal elements of the result can be obtained 
simply by multiplying the two corresponding diagonal elements. This results that the trace of the £ on a periodic 
orbit is the following: 



tr Cp - tr Ls^ is2 • • • is„p - 2^ I A Um+i) 

m=0 I^^Pl 

_ 1 ^ 1 _ 1 _ 1 

~ M,^„ AT " |A,| (1- 1/A,) - ^^A;| (^^^ 

The absolute value is a consequence of the sign(/') factor in (pO|). For the total trace, using ( pi] ) we get back the 
standard deterministic trace formula [Q 

tr/:" = ^np^,5„,„^rT— — TTj (36) 

p r=l I P' 



C. Numerical tests 

Here we continue the calculations of Sect. ? of ref. |g], where more details and discussion may be found. We test 
our perturbative expansion on the repeller of the l-dimcnsional map 

^1 /I ^*' 



/(.) = 20(^--^--.Jj. (37) 

This repeller is a nice example of an "Axiom A" expanding system of bounded nonlinearity and complete binary 
symbolic dynamics, for which the deterministic evolution operator eigenvalues converge super-exponentially with the 
cycle length. 

The inverse of / has two branches: 



r\y) = 



2 lie 20/ '2 lie 20 



(38) 



Due to the symmetry of (^^ around x = 0.5, the derivatives on the two inverse branches have the same absolute 
value, but opposite signs. We compute the leading eigenvalue of the evolution operator (the repeller escape rate) in 
the presence of Gaussian noise, using three complementary approaches. The perturbative result in terms of periodic 
orbits and the weak noise corrections is compared to the eigenvalue computed by a numerical lattice discretization in 
ref. |l|] . In the preceding paper (ref. ) we compared the numerical eigenvalue with the a^ result and estimated the 
coefRcient of a^ to be approximately 2700. Here we compute the order a^ coefficient 2076.47 ... to 14 digits accuracy, 
as well as the u* and ct^" coefficients. Furthermore, we estimate the asymptotic form of the ct^™ coefficient. 
The numerical calculations of tr£" proceeds as follows: 

1. Creating periodic orbits up to length 10 using the method of iterating backwards Q 

2. Computing [24 x 24] B matrices at each point of orbits 

3. From the B matrices computing [16 x 16] L matrices at each point of orbits, up to cr^" corrections in each 
matrix-element 

4. Multiplying the L matrices. At shorter orbits repeating this r times. (rrip — n) 

5. Computing the traces of the result matrices, multiplying them by Up and adding them up to get tr£" as in ( |ll| ) 

6. Create cummulants from the traces. 

7. Find zeros. 

At n = 2 [20 X 20] sized L (and corresponding to that [28 x 28] sized B) matrices were used because of slower 
convergence in tr£" while enlarging the matrix size. At n = 1 [26 x 26] is the size in L one has to reach to get decent 
result for tr C. 



D. From traces to cr "* corrections 

To get the leading eigenvalue of the Perron- Frobenius operator from the traces of £", we follow the method outlined 
in y]. The trace of £" for n — 1,2, .... 10 can be written as: 

oo 

tr/:"=^C„Ja^ (39) 

3=0 

where we know the value of the first five non-zero C„o, ■ • • , C'n.io ■ The cumulants Qnj in 

oc oo 

det(l -zC) = l-^Y^ Q„j2"crJ' (40) 

are obtained recursively as 

_. / n—l m \ 

Qnm = ~ I ^mn ~ / ^ / ^ Qk,7n-l'-^n — k,l I • (41) 

"" \ fe=l (=0 / 

For Gaussian noise only even / terms contribute. Let zq be the solution of the noiseless condition det(l — zC)\^^q = 0. 
Expanding the spectral determinant around zq and a^ — 0, we write: 

det(l - (z + zo)C) = F- Fioz - Fo2<J^ - F2oz^ - Fi2zcr^ - Fo^a'^ 

-0{z^,(j'^) (42) 

where the coefficients can be obtained from the cumulants as: 

F = 1 — J2m=l QmfiZ^ Fio — J2m.= l ^Q-mfiZ"^ 

F02 = El^lQn.a^S' F20 = iE™=2"l("^-l)Q™.0^0""' (43) 

F12 — J2rn=l ^^Qrn.2Z^ Fo4 = J2m=l Qrn,4Z^ 

If we want to calculate only the a corrections just up to fourth order, first we have to solve this equation: 

F - Fio(z2cr2 + z^a^) - Fo2(t2 - ^20(^2^^ + z^a^f 

-F12 {z2a^ + Z4CT*)ct2 - Fo4 cr-* = 0. (44) 

Since cr is a parameter, the coefficients of ct*', or a^ or a* have to vanish: 

F = 

—F10Z2 — F02 = 
^2 



'F10Z4 — F2aZ2 — F12Z2 — Foi — 



The first equation gives solution for zq, which numerically was done by Newton's method. From the second we can 
express Z2, from the third Z4: 

-F02 F2QZ2 + F12Z2 + Fq4 

22 = -T^ Z4 = ^ (45) 

-rio -rio 

To get the eigenvalue we have to reciprocate z: 

2 , 4 1 1 2:2 2 fz4 ■' 

vq + V20 + v^a = ; 5— = 2 cr - -T - - 

ZO + 220- + Z4Cr* Zq Z^ \ Zq ; 

+0(a*^) (46) 



We see that: 



1 22 Zi zl 

Vq = —, V2 = 2' ^4 = -— + 73' (47) 



3 ■ 

Zo ^0 ^0 -^0 



or by substituing into these mW: 



vo 



1 



^f 



V2 



02,, 2 



F^ 



10 



0' 



i^A — -^{F2qFq2 — F12F02F10 + Fo^Fif^ + i^QjFio 1^0)^0. 



(48) 



To go up further in the powers of a corrections, we have to do the same as we have done for order four, the details of 
this is in the appendix. 



1^0 



V2 



Vi 



0.308 

0.37140 

0.3711096 

0.371110995255 

0.371110995234863 

0.371110995234863 



0.42 

1.422 

1.43555 

1.435811262 

1.43581124819737 

1.43581124819749 



2.2 
32.97 
36.326 
36.3583777 
36.35837123374 
36.358371233836 



I'd 



i^s 



17.4 

1573.3 

2072.9 

2076.479 

2076.4770492 

2076.47704933320 

2076.47704933321 



168.0 

112699.9 

189029.0 

189298.8 

189298.12802 

189298.128042526 

189298.128042526 



TABLE I. Significant digits of the leading deterministic eigenvalue and its a 



, a coefficients, calculated from the spectral 



determinant as a function of the cycle truncation length n. Note the super-exponential convergence of all coefficients. We have 
computed all cycles to length 10, but contributions of those longer than n = 6 lie below the machine precision. 



The perturbative corrections to the leading eigenvalue (escape rate) of the weak- noise evolution operator are given in 
table H, showing super-exponential convergence with the truncation cycle length n. The super-exponential convergence 
has been proven for the deterministic, vq part of the eigenvalue |§,0, but has not been studied for noisy kernels. It 
is seen that a good first approximation is obtained already at n — 2, using only 3 prime cycles, and n = 6 (23 prime 
cycles in all) is in this example sufficient to exhaust the limits of double precision arithmetic. The exact value of 
vq = 2076.47 ... is not wildly different to our previous numerical estimate [^ of 2700. 

VI. TESTING THE RESULTS 

In previous papers the eigenvalue was computed numerically. In fig. 0to fig. we check how well our results fit that. 
The present result fits the data of the numerical discretisation method well in the small noise region, but not above 
sigma = 0.1. At given sigma, a matrix representation of the Perron-Frobenius operator can be obtained numerically, 
at each matrix element carrying out a numerical integration: 



L 



l.k 



d' 
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I dze- 


(b-/(^))^ 


dyi 


.\/2^afc! J 





(49) 



y=0 



The leading eigenvalue of this matrix fits well the discretisation data well at large sigma, as it is shown on fig. 

The cumulants Q^asafunctionsofcycle length n 
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FIG. 1. The generalized cumulants Q„j as a function of cycle length n for j = 0,2,4,6,8. Qno is the cumulant of the 
noiseless case. Superexponential convergence can be observed until cycle length n — 6, then numerical errors take over. 
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The numerical eigenvalue and known terms 
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FIG. 2. The numerical eigenvalue and the known ternis:j^(n) = X/fc=o ^'^^cr 
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The difference between numerical eigenvalue and known tenns 
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FIG. 3. The difference between the numerical eigenvalue and the known terms. 
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FIG. 4. The result of numerical integration and the eigenvalue fit together well in the large a region. 



VII. SUMMARY AND OUTLOOK 

We have formulated a perturbation theory of stochastic trace formulas based on a polynomial basis matrix repre- 
sentations, expanded around infinitely many chaotic saddle points (unstable periodic orbits). 

We note in passing that for \-d repellers a diagonalization of an explicit truncated Cmn matrix yields many more 
eigenvalues than the cycle expansions |Iq,10| . The reasons why one persists anyway in using the periodic orbit theory 
are partially aestethic, and partially pragmatic. Explicit Cmn demands explicit choice of a basis and is thus non- 
invariant, in contrast to cycle expansions which utilize only the invariant information about the flow. In addition, we 
do not know how to construct £,„„ for a realistic flow, such as the 3-disk problem, while the periodic orbit formulas 
are general and straightforward to apply. 
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APPENDIX A: ALGEBRA 

To go up further in the powers of a corrections in i/, we have to do the same as we have done for order four. First 
expand (|4|): 

det(l - (z + zo)C) =F- Fwz - Fmcr'^ - F2az^ - F^za^ - Fo^a"^ 
-Fsoz^ - F22Z^cr^ - FuZG^ - Foea^ 



-Fioz^ - FaazV^ - F24ZV* - Fi^za"^ - Fosa"^ 
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The coefficients are: 
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The equation to be solved is: 



-Fi2{z2(T^ + 2;4ct'' + Zgcr^ + Z^(T^)(T^ - Fo4Cr'' - ^30(22^^ + Z4Cr'' + ZgCr^ + Z^a^f 
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The solutions for zg and zs are: 
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The connection between the v^s and the z's: 
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We can say that: 
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